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Surface  heat  transfer  and  boundary  layer  transition j measurements 
have  been  performed  on  a  slender  blunted  cone  “  8.0. 

These  measurements  have  been  obtained  on  the  conrcal  surface  with 
mass  transfer  occurring  on  the  spherical  nose  cap.  The  free 
stceaniyReynold*-nwrtber ,  based  on  nose  radius,  varied  between 
l  p  .  6x10*/  and  <5\  9x10 zJ.  Data  was  obtained  both  at  zero  angle  of  at¬ 
tack  and  at  angles  )as  high  as  14 f ,  which  exceeds  the  semi-vortex  J 
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cone  angle  of  10*.  'The  technique  utilized  in  the  determination 
of  transition  location  was  a  thin  film  gauge,  on  the  model  surface, 
which  was  maintained  at  a  constant  temperature  by  an  anemometer 
system.  The  RMS  output  of  the  film  gauge  was  found  to  accurately 
locate  the  onset  of  boundary  layer  transition. 

A  numerical  analysis  of  film  cooling  in  the  stagnation  region  of  a 
slender  cone  was  also  performed  using  the  thin  shock  layer  approx¬ 
imation.  The  porous  region  was  analyzed  and  compared  to  pre¬ 
viously  obtained  results  using  different  assumptions.  The  down¬ 
stream  region,  at  angle  of  attack,  utilize  the  porous  region  re¬ 
sults  as  initial  conditions.  The  computer  program,  however,  was 
not  successfully  operated  in  this  mode. 
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POLY  M/AE  Report  No.  82-10 


Summary 

The  studies  performed  under  the  sponsorship  of  this  grant 
are  both  of  an  experimental  and  an  analytical  nature;  those 
are  summarized  below: 

(i)  The  experimental  portion  of  tho  study  involved  the 

measurement  of  transition  initiation  on  slender  blunted 

cones  in  a  hypersonic  flow  with  mass  transfer  cooling. 

This  was  done  for  both  zero  and  non- zero  angle  of  attack 

conditions  and  included  the  development  of  a  transition 

detection  technique  which  utilized  a  thin  surface  film 

gauge.  A  correlation  of  Re„  was  obtained  which  is 

ytr 

useful  in  the  prediction  of  transition  under  thpBe  con¬ 
ditions.  The  other  aspect  of  the  experimental  program 
was  the  utilization  of  laser  diagnostics  to  measure  li¬ 
quid  water  concentrations  in  a  film  cooled  boundary 
layer.  These  programs  are  discussed  in  more  detail  in 
the  following  section. 

(ii)  The  theoretical  portion  of  the  program  centered  on 
the  numerical  analysis  of  film  cooling  in  the  stagnation 
region  of  a  slender  cone,  utilizing  the  thin  shock  layer 
approximation.  The  porous  region  was  completely  an- 
alzed  and  compared  to  previous  analysis  using  different 
assumptions;  this  resulted  in  a  technical  publication. 

The  downstream  region  was  also  analyzed  and  successfully 
tested  at  zero  angle  of  attack.  The  angle  of  attack  re¬ 
sults  would  utilize  the  porous  region  results  as  Initial 
conditions  and  are  also  expected  to  produce  useful  infor- 
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mation  although  the  computor  program  was  not  tested  in 
this  mode. 

The  principal  investigators  on  this  grant  for  the  period 
October  i.980  through  March  1982  were  Prof.  M,  H.  Bloom  and 
Prof.  R.  J.  Cresci. 
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TECHNICAL  DISCUSSION 

On  slender  bodies  traveling  at  hypersonic  speeds,  whether 

mass  is  injected  by  an  active  coolant  system  or  whether  the 
injection  occurs  as  a  result  of  surface  ablation,  the  major 
transfer  of  mass  will  occur  in  the  stagnation  region  followed 
by  a  surface  on  which  the  effective  injection  rate  is  effec¬ 
tively  zero.  In  the  ablation  system,  this  is  caused  by  the 
radiation  from  the  high  temperature  gases  behind  the  bow  shock 
wave  in  the  nose-cap  region.  As  the  shock  weakens  in  the  down¬ 
stream  direction  and  becomes  more  conical  in  shape,  the  gas 
temperature  drops  rapidly  and,  along  with  it,  the  attendant  sur¬ 
face  ablation.  If  one  considers  the  active  coolant  system, 
again  one  finds  that  the  surface  mass  transfer  is  confined  to  the 
stagnation  region  since  in  this  case  the  minimum  injection  re¬ 
quired  for  reduction  of  conveotive  heating  is  the  important  con¬ 
sideration.  In  either  situation,  however,  there  will  be  some 
residual  cooling  effect  as  the  injected  fluid  flows  downstream 
within  the  boundary  layer. 

It  has  been  observed  in  ref.  (1)  that  this  increased  boundary 
layer  mass  flow  will  initiate  transition  to  turbulence  earlier 
than  in  a  boundary  layer  with  no  upBtream  injection.  In  a  trans¬ 
itional  or  turbulent  boundary  layer,  the  local  heat  transfer  is 
usually  significantly  higher  than  that  in  a  laminar  flow  at  the 
same  free  stream  conditions.  As  a  result,  it  is  extremely  impor¬ 
tant  to  be  able  to  determine  the  onset  of  transition  when  such  ac¬ 
tive  cooling  systems  are  proposed,  or  when  they  occur  naturally 
through  surface  ablation.  Downstream  effects  of  mass  injection  on 
a  slender  cono  were  obtained  in  reference  (2)  from  surface  heat 
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transfer  measurements?  these  measurements  were  then  used  to  attempt 
to  correlate  the  transition  location  with  mass  transfer  rate  and 


free  stream  Reynolds  numbers.  This  iB  the  only  known  previous 

*  1 

study  of  transition  located  on  the  impermeable  surface  downstream 
of  the  mass  transfer  region. 

The  present  study  deals  with  a  combined  experimental  and  theo¬ 
retical  investigation  of  the  downstream  effects  of  mass  transfer 
in  the  stagnation  region  of  slender,  conical  bodies,  of  particular 
interest  is  the  transitional  behavior  of  the  boundary  layer  in  the 
presence  of  upstream  mass  injection.  Both  zero  and  non-zero  angles 
of  attack  were  considered  in  both  the  experimental  and  theoretical 
portions  of  the  study.  Reference  (3)  describes  both  portions  of 
the  effort  in  a  reasonably  detailed  fashion. 

Experimental  studies 

Ons  purpose  of  the  present  study  was  to  attempt  to  measure  the 
onset  of  transition,  on  the  impermeable  surface,  with  greater  ac¬ 
curacy  than  is  possible  by  interpretation  of  surface  heat  transfer 
measurements.  This  was  accomplished  by  using  thin  surface  films 
which  were  heated  and  monitored  by  a  constant  temperature  hot  wire 
anemometer  system.  This  technique  was  used  to  accurately  locate 
the  transition  point  under  a  variety  of  free  stream  Reynolds  num¬ 
bers,  surface  locations,  and  mass  injection  rates.  The  study  was 
performed  at  a  free  stream  Mach  number  of  8.0  on  a  spherically 
blunted  oone  of  10°  half  angle,  as  seen  in  figure  (1);  reference 
(4)  describes  this  study  in  detail.  Although  the  major  portion  of 
this  study  was  conducted  at  zero  angle  of  attack,  some  tests  were 
run  at  angles  of  attack  between  2 6  and  14°.  Both  the  downstream 


of f ectivenoss  and  the  transition  location  wore  also  measured  undor 
these  flow  conditions.  The  testa  were  performed  in  a  Mach  8.0 
blow  down  wind  tunnel,  which  has  a  two  foot  diameter,  axisymme- 
tric  test  section.  The  freestream  stagnation  temperature  of  the 
flow  was  maintained  at  2000°R.j  this  is  accomplished  by  heating 
the  air  in  a  pebble  bed  heater.  The  stagnation  pressure  of  the 
air  in  the  heater  may  be  varied  from  50  to  600  psia.  The  above 
stagnation  conditions  correspond  to  a  range  of  freestream  Reynolds 
numbers  of  1.13xl04  to  1.35x10s  per  inch. 

Although  surface  heat  transfer  measurements  were  obtained  in 
the  current  test  program,  the  primary  objective  was  to  develop  a 
technique  which  more  accurately  ascertained  the  transition  loca¬ 
tion  on  the  model.  In  the  past,  various  methods  have  been  uti¬ 
lized  in  the  determination  of  transition  loaation,  e.g.  ,  surface 
heat  transfer,  Preston  tubas,  and  thin  film  gauges.  In  the  cur¬ 
rent  test  program  it  was  decided  to  use  the  thin  surface  film 
technique,  which  utilizes  standard  TSI  miniature  flush  mounted 
sensors  in  which  the  film  is  deposited  on  the  end  of  a  quartz 
cylinder  1.5mm  in  diameter.  The  gauge  is  maintained  at  a  tempera¬ 
ture  slightly  higher  than  that  of  the  surrounding  model  by  a  con¬ 
stant  temperature  anemometer  which  is  used  to  power  the  gauges 
and  obtain  the  RMS  sensor  output. 

The  output  of  the  anemometer,  the  RMS  meter,* the  thermo¬ 
couples  and  the  pressure  transducers  was  fed  into  a  bank  of  Honey¬ 
well  Aooudata  122  linear  amplifiers.  The  amplifiers  normalized 
all  incoming  signals  to  a  range  of  0.0  to  1.0  volts.  This  normal¬ 
ized  output  was  reoordod  on  a  Honeywell  Model  101  AM-FM  magnetic 
tape  recorder.  The  microprocessor  of  the  Model  101  can  mark  each 
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channel  with  a  timing  pulse  so  that  all  recorded  data  may  be  cor¬ 
related  In  real  time. 

After  the  wind  tunnel  teat  run  is  terminated,  the  data  re¬ 
corded  on  the  tape  recorder  is  played  back  into  a  Digital 
PDP11-34  computer  which  has  an  AR-11  laboratory  peripheral 
system.  The  AR-11  system  is  capable  of  analog  to  digital  conver¬ 
sion  of  16  channels  of  analog  data,  which  is  stored  by  the 
PDP11-34  on  a  diskette  and  is  then  ready  for  processing  and  anal¬ 
ysis. 

In  order  to  obtain  transition  identification  in  a  single 
test  run,  the  free  stream  test  Reynolds  number  was  held  constant 
while  the  masB  transfer  rate  was  varied.  The  results  of  a  typi¬ 
cal  test  are  shown  in  figure  (2) ,  which  presents  the  coolant 

chamber  pressure  vs.  time  and  the  RMS  signal  vs.  time  at  a  surface 

—  4 

location  of  s  -  36.8  and  a  free  stream  Reynolds  number  of  4.0x10  . 

The  location  of  transition  is  clearly  evident  on  the  RMS  plot  and 
determines  the  injected  mass  transfer  rate  required  to  initiate  ' 
transition  at  the  specified  surface  location  and  at  the  test  Rey¬ 
nolds  number  chosen. 

These  data  were  correlated  for  a  wide  range  of  surface  loca¬ 
tions,  Reynolds  numbers  and  mass  transfer  rates.  The  results  of 
this  correlation  are  shown  in  figure  (3)  in  which  all  the  data 
are  seen  to  follow  the  same  trend  in  terms  of  the  *  transitional 
Reynolds  number  based  on  momentum  thickness. 

One  of  the  more  interesting  results  of  these  tests  was  that 
transitional  flow  was  observed  close  to  the  injection  region,  fol¬ 
lowed  by  a  region  of  laminar  flow  which  in  turn  was  followed  by 

< 

another  region  of  transitional  and  ultimately  turbulent  boundary 


layer.  This;  behavior  occurred  only  at  certain  Reynolds  numbers 
and  implied  a  local  rolaminarization  of  the  boundary  layer. 

Figure  (.4)  presents  the  mass  transfer  required  to  obtain  transi¬ 
tion  as  a  function  of  surfaoe  location  and  free  stream  Reynolds 
number.  The  data  obtained  at  the  lower  Reynolds  number  is  what 
is  typically  expected,  however,  at  ReM  »  0.4xl05  one  may  observe 
a  transition  reversal  that  occurs  at  values  of  7  less  than  25. 

At  angle  of  attack,  the  surface  film  gauge  UBed  to  detect 
boundary  layer  transition  was  only  run  at  the  I  -  9.09  location. 
This  data  is  shown  in  figure  (5)  for  three  meridian  planes t 
windward  (.$  ■  0)  cross  plane  (4>  ■  90°)  ,  and  one  intermediate 
plane  (i  -  60°). 

Some  general  conclusions  can  be  reached  from  this  test  pro¬ 
gram  and  may  be  summarized  as  follows i 

(i)  An  accurate  determination  of  the  transition  loca¬ 
tion  may  be  obtained  by  using  thin  surface  film 
gauges,  operated  by  a  constant  temperature  anemometer 
Bystem,  by  examination  of  the  RMS  signal, 
tii)  A  correlation  of  transition  Reynolds  numbers,  based 
on  momentum  thickness,  has  been  obtained  in  terms  of 
the  surface  location  and  the  mass  transfer  required 
to  initiate  transition. 

'  (.iiil  For  a  certain  range  of  free  stream  Reynolds  numbers, 
the  effect  of  mass  transfer  was  to  initiate  early 
transition,  followed  by  a  region  in  which  relaminar- 
ization  of  the  boundary  layer  occurred. 

Civ)  Both  surfaoe  heat  transfer  'and  transition  data  were 

obtained  up  to  angles  of  attack  greater  than  the  half 


cone  angle.  Reasonable  mass  injection  rates  were 
found  to  effectively  reduce  the  nose-tip  heat  trans¬ 
fer  to  values  which  exist  on  the  downstream  conical 
surface,  even  at  the  highest  angles  of  attack  con¬ 
sidered. 

(v)  On  the  windward  meridian,  the  presence  of  non-zero 
angle  of  attack  delayed  boundary  layer  transition, 
even  in  the  presence  of  the  destabilizing  mass 
transfer. 

Theoretical  Studies 

Reference  (5)  presents  a  special  case  of  the  more  generalized 
equations  developed  in  reference  (3) .  For  this  analysis  a  limiting 
form  of  the  full  axisymmetric  shock  layer  equations  is  solved  in 
the  stagnation  region.  The  limiting  form  of  the  equations  ignores 
the  elliptic  nature  of  the  more  general  equations  and  is  equivalent 
to  assuming  local  similarity.  The  streamwise  momentum  and  conti¬ 
nuity  equations  are  evaluated  using  the  coupled  block  tridiagonal 
algorithm.  Energy  is  solved  uncoupled  from  the  rest  of  the  system* 
the  influence  of  temperature  is  less  critical,  as  the  characteris¬ 
tic  shape  of  the  velocity  profiles  is  governed  by  the  streamwise 
momentum  and  continuity  equations.  A  variable  grid  is  employed 
that  ensures  maximum  resolution  about  the  viscous  interface. 

Solutions  to  the  system  of  equations  have  been  found  for  a 
range  of  free-stream  Mach  numbers,  Reynolds  numbers,  and  injection 
rates  up  to  250%  of  the  free-stream  mass  flux  per  unit  area.  The 
largest  injection  rates  available  in  the  literature  appear  to  ap¬ 
proach  80%  of  the  free-stream  mass  flux.  Zero-injection  solutions 
for  the  present  analysis  are  also  presented  in  order  to  show  the 
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range  and  consistency  of  the  approximations.  At  low  and  moderate 
ihjection  rates,  comparisons  with  experimental  data  and  alterna¬ 
tive  theories  are  given  for  profiles,  standoff  distance,  and  in- 

i  * 

terface  location.  Effects  of  grid  size  are  discussed.  Compari¬ 
sons  of  stagnation  region  profiles  with  downstream  profiles  are 
discussed,  as  is  utilisation  of  the  solution  as  initial  conditions 
for  more  complete  numerical  analyses,  which  could  determine  the 
downstream  effects  of  large  localized  upstream  injection. 

The  following  conclusions  may  be  obtained  regarding  this 
studyi  at  low  injection  rates,  solutions  compare  well  with  a 
variety  of  previous  theories  and  experimental  data.  For  very 
large  injection  rates  the  statio  pressure  variation  is  significant 
and  must  be  taken  into  account.  Even  at  high  injection  rates  and 
large  Reynolds  numbers  good  resolution  can  be  achieved  with  a 
moderate  number  of  points  if  the  equations  are  dast  in  conserva¬ 
tion  form  and  streamwise  momentum  and  continuity  are  solved  as  a 
ooupled  system.  The  number  of  iterations  necessary  for  conver¬ 
gence  varied  from  12  for  the  zero- injection  aase  at  low  Reynolds 
numbers  to  35  for  the  injection  case  at  high  Reynolds  numbers. 

This  is  not  excessive  considering  the  very  large  gradients  in¬ 
volved  in  the  viscous  interface. 

The  solutions  are  ideally  suited  as  initial  conditions  for 
marahing  teohniques.  Any  downstream  influence  on  the  stagnation 
solution  may  be  removed  by  global  iteration.  The  algorithm  can 
be  easily  applied,  is  stable,  and  handles  large  gradients  without 
difficulty.  Downstream  profiles  compare  well  with  the  stagnation 

solution  in  nondimensional  shock  layer  coordinates.  Keeping  this 

» 

in  mind,  it  might  be  possible  to  initialize  a  time-dependent  tech- 
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nique,  which  requires  initial  values  for  the  entire  flow  field, 
knowing  only  the  stagnation  region  solution. 

This  analysis  was  developed  to  generate  initial  conditions 
for  a  more  complete  three-dimensional  analysis  of  the  downstream 
effect  of  large  localized  upstream  injection  on  hypersonic  blunt 
bodies . 

The  downstream  region  of  interest  in  the  overall  program 
was  also  treated  by  a  numerical  analysis  of  the  shock  layer  equa¬ 
tions.  This  analysis  was  set  up  to  include  both  zero  and  non¬ 
zero  angles. of  attack.  The  principal  purpose  of  the  analysis 
was  to  develop  a  numerical  solution  to  the  shook  layer  equations 
reflecting  the  effects  of  large  upstream  injection.  To  deal  with 
the  severe  gradients  characteristic  of  large  mass  injection, 
numerical  techniques  are  utilized  which  have  been  shown  to  be 
capable  of  computing  difficult  flow  problems.  .  An  unconditionally 
stable  blook-tridiagonal  algorithm  provides  a  coupled  solution  of 
the  continuity,  streamwise  momentum,  and  transverse  momentum  equa¬ 
tions.  A  modified  central  space  difference  equation  is  used  for 
all  convective  transverse  first  derivatives  of  the  transverse 
velocity  to  insure  the  unconditional  stability  of  the  algorithm. 

The  energy  equation  is  Bolved  uncoupled  from  the  rest  of  the 
■yatem  by  a  two-dimensional  strongly  implicit  algorithm,  which  is 
capable  of  updating  the  solution  to  the  energy  equation  in  the  en¬ 
tire  transverse  plane  in  one  sweep.  These  techniques  require  some¬ 
what  more  computer  storage,  however,  they  provide  the  necessary 
coupling  of  the  boundary  conditions  and  the  dependent  variables 
for  a  high  degree  of  efficiency  and  rapid  convergence  rates. 

To  provide  the  upstream  mass  injection  a  porous  injection 


region  subtending  a  spherical  solid  angle  of  60°  is  loaatod  at 
the  apox  of  the  sphoro-cone  configuration  (fig.  1) .  At  present 
no  provision  for  cross  flow  separation  is  incorporated  into  the 
analysis.  Initial  conditions  are  provided  for  angle  of  attack 
oases,  by  solving  tho  fully  viscous  shock  layer  equations  at  an 
axisyrametric  stagnation  point  for. large  rates  of  injection.  The 
stagnation  solution  is  then  rotated  to  a  body  fixed  coordinate 
system  to  initialise  the  three  dimensional  solution.  The  method 
requires  no  initial  guess  for  the  entire  shook  shape.  The  only 
assumption  made  on  the  shook  is  when  the  stagnation  region  ini¬ 
tial  profiles  are  computed,  otherwise  the  shook  shape  is  pre¬ 
dicted  at  each  step  in  the  computation.  The  Rankine-Hugoniot 
conditions  are  used  to  compute  the  outer  boundary  conditions  on 
the  shock  layer.  Tha  sphere  cone  interface  is  treated  as  a  dis¬ 
continuity  in  curvature,  therefore,  no  differences  in  the  flow 
variables  are  taken  across  the  sphere-cone  juncture.  Provision 
is  made  in  all  the  difference  equations  for  variable  grids  in 
all  coordinate  direotions.  This  is  important  in  thiB  case  for 
the  mass  injection  is  abruptly  terminated  at  a  particular  down¬ 
stream  location.  The  streamwise  grid  spacing  at  the  injection 
out  off  point  must  be  fine  enough  to  handle  the  streamwise 
transition  from  injection  profiles  to  non-injection  profiles. 

At  this  location  discontinuities  in  surface  heat  transfer  and 
shear  stress  have  been  found  to  occur. 

It  is  realised  that  choosing  the  three  dimensional  coupled 

algorithm  to  solve  the  two  momentum  equations  and  the  continuity 

» 

equation,  as  well  as  choosing  the  strongly  implicit  plane  algori- 

*  * 

thm  to  solve  the  energy  equation  requires  more  computer  storage 


9 


than  a  local  cascading  scheme.  This  increase  in  computer  storage 
is  outweighed  by  the  stability  of  the  methods,  the  reduction  in 
the  number  of  iterations,  the  amount  of  grid  points  necessary 
and  a  possible  increase  in  the  accuracy  of  the  solution. 

The  details  of  this  analysis  are  described  more  fully  in  ref¬ 
erence  (3),  however,  no  numerical,  results  were  available  for  the 
angle  of  attack  condition.  Computer  operating  difficulties  and 
the  posslblity  of  undiscovered  bugs  in  the  program  prevented  the 


successful  operation  of  ths  program  prior  to  the  expiration  of 
the  grant. 

Laser  Diagnostics 

Ths  development  of  an  optical,  nonintrusive  flow  diagnostic 
system  using  laser  Raman  scattering  techniques  has  also  been  un¬ 
derway  to  determine  water  vapor  and  liquid  droplet  concentrations 
in  a  high  spaad  boundary  layer.  The  presenoe  of  the  water  orig¬ 
inates  whan  a  liquid  film  is  utilized  for  nose  tip  cooling.  As 
the  film  progresses  downstream,  it  breaks  up  into  droplets  (or 
vaporizes)  and  it  is  therefore  important  to  know  the  particular 
state  of  the  water  at  any  location  in  order  to  be  able  to  evaluate 
the  physical  mechanism  of  film  deterioration.  In  its  original 
form,  the  system  waB  based  on  the  continuous  wave,  low  level,  co¬ 
herent  anti-Stokes  Raman  Scattering  (CARS)  arrangement.  This 
system  employing  two  argon-ion  lasers  had  to  be  abandoned  due  to 
the  nonreliability  of  the  lasers  supplied  to  us  by  the  Control 
Laser  Corp.  of  Orlando,  Florida.  A  new  laser  system  whioh  became 
available  to  us  recently  has  therefore  been  utilized  for  the  de¬ 
velopment  of  a  reliable  CARS  system  to  measure  water  vapour  oon- 

*  * 

centrations.  Results  have  been  obtained  by  this  sytem  whioh  is 
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capable  of  resolving  relatively  minor  concentrations  of  water 

vapour. 

An  effort  was  made  to  obtain  the  same  results  for  liquid  water 
as  for  the  vapour  in  the  laboratory  using  the  same  laser  system. 
Normally,  this  would  require  either  additional  tests  or  another 
laser  system,  however,  an  attempt,  was  made  to  obtain  the  concen¬ 
trations  of  liquid  water  and  water  vapour  simultaneously  using 
only  one  laser.  This  system  would  allow  one  to  determine  the  com¬ 
plete  phase  balance  of  the  coolant  at  any  instant  of  time,  or  at 
any  location  along  the  surface. 

As  part  of  the  development  of  the  laser  diagnostic  system  for 
the  simultaneous  acquisition  of  concentration  of  water  vapour  and 
liquid  water  droplets,  it  has  been  possible  to  obtain  both  with 
the  same  laser,  using  separate  CARS  systems.  This  process  utiliz¬ 
ing  the  single.  Neodynium  Yag  laser  and  the  appropriate  specially 
designed  tunable  dye  laser,  derived  from  and  pumped  toy  a  fraction 
of  the  main  laeer  power,  has  been  halted  due  to  the  lack  of  time. 

It  is,  however,  believed  that  the  approach  is  fundamentally  sound 
and  should  be  pursued  in  the  future. 
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FIG  3  CORRELATION  OF  TRANSITIONAL  Ren  WITH  MASS  TRANSFER 
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FIG  5  TRANSITIONAL  BEHAVIOR  at  $=9.09  vs. 
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